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Abstract 

Determining the spin and the parity quantum numbers of the recently 
discovered Higgs-like boson at the LHC is a matter of great importance. In 
this paper, we consider the possibility of using the kinematics of the tagging 
jets in Higgs production via the vector boson fusion (VBF) process to test 
the tensor structure of the Higgs-vector boson (HVV) interaction and to 
determine the spin and CP properties of the observed resonance. We show 
that an anomalous HVV vertex, in particular its explicit momentum depen- 
dence, drastically affects the rapidity between the two scattered quarks and 
their transverse momenta and, hence, the acceptance of the kinematical cuts 
that allow to select the VBF topology. The sensitivity of these observables 
to different spin-parity assignments, including the dependence on the LHC 
center of mass energy, are evaluated. In addition, we show that in associ- 
ated Higgs production with a vector boson some kinematical variables, such 
as the invariant mass of the system and the transverse momenta of the two 
bosons and their separation in rapidity, are also sensitive to the spin-parity 
assignments of the Higgs-like boson. 



1. Introduction 

The mechanism for electroweak symmetry breaking ll[ is minimally re- 
alized in the Standard Model (SM) of particle physics [2( by introducing a 
doublet of complex scalar fields that develops a non-zero vacuum expectation 
value; this leads to the existence of only one physical scalar boson H, with 
the spin, parity and charge conjugation assignments of a J PC = ++ state. 
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Accumulated data by the ATLAS and CMS collaborations seem to indi- 
cate that the newly observed bosonic particle with a mass around 125 GeV 
[3| has production and decay rates that are compatible with those expected 
for the SM Higgs particle. However, it will be necessary that experiments 
test the spin and CP properties of the new boson before we can truly identify 
the observed particle with the Higgs boson of the minimal SM. In particular, 
more data will be required in order to verify that the Higgs sector is not 
extended in order to contain CP-violating interactions that could make the 
observed Higgs particle a mixture of CP-even and CP-odd states. In addi- 
tion, while it is clear from the observation of the 77 decay channel that the 
observed bosonic state cannot have spin equal to one j3| , it is not yet entirely 
excluded that it has a spin equal to two or more. 

It is well known that the structure of the Higgs to vector boson (HVV) 
coupling is a possible tool to probe the Higgs J PC quantum numbers at the 
LHC and the ILC. Distributions in various kinematic variables, the virtuality 
of the off-shell gauge boson, and angular correlations of the decay products 
in the decay channel H — > VV — > 4f with V = W, Z in particular, allow to 
discriminate between CP-even and CP-odd states as well as spin-zero from 
higher spin particles [El, 0, 0, @, @] • It is also known that the azimuthal angle 
distribution of the two outgoing forward tagging jets in Higgs production in 
the vector boson fusion (VBF) process effectively discriminates between the 



CP quantum numbers of a scalar resonance [10|, III • 

In the present paper, we consider Higgs production in the vector boson 
fusion mechanism, pp — > Hjj, in the presence of an anomalous Higgs-vector 
boson vertex that parametrises different spin and CP assignments of the 
produced state. The anomalous HVV coupling is introduced by allowing for 
an effective Lagrangian with higher dimensional operators, that include four 
momentum terms which are absent in the SM. We show that the kinematics 
of the forward tagging jets in this process is highly sensitive to the structure 
of the anomalous HVV coupling and that it can effectively discriminate 
between different assignments for the spin (spin-0 versus spin-2) and the 
parity (CP-even versus CP-odd) of the produced particle. 

We find, in particular, that the correlation between the separation in 
rapidity and the transverse momenta of the scattered quarks in the VBF 
process, in addition to the already discussed distribution of the azimuthal 
jet separation, can be significantly altered compared to the SM expectation. 
These kinematical variables define new corners of the phase-space that have 
not been explored by the experiments at the LHC to probe anomalous HVV 
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couplings and check the J PC assignments of the newly observed particle. 
Some of these observables significantly depend on the center of mass energy 
and strong constraints on anomalous couplings can be obtained by perform- 
ing measurements at the LHC with energies of y/s = 8 TeV and y/s = 13 
TeV. Finally, we also consider associated Higgs production with a massive 
gauge boson , qq — > HV (VH, V = W ± and Z ), and show that the invariant 
mass of the VH system as well as the transverse momenta and rapidities of 
the H and V bosons are also sensitive to anomalous HVV couplings. 

An analysis that is, in some aspects, similar to ours has been recently 



reported in Ref. [12[ and we will discuss the main differences between this 
study and the one presented here at the end of the paper. 

For the rest of this article, the next section describes the physical set-up 
and sections [3] and H] our phenomenological analyses in the VBF and Higgs- 
strahlung processes; section [5] summarizes the main findings. 

2. The physical set-up 

In the SM, the couplings of the Higgs boson to the massive electroweak 
gauge bosons are precisely formulated and come out as gHvv ^ gM v V^V^ 
with g the SU(2) coupling constant. However, this is not the most gen- 
eral form of the Higgs-gauge boson vertex. Parametrising the coupling of 
a scalar state to two vector bosons in the form iT^ u (p, q)€ fl (p)el(q), one 
can write down the most general form of the HVV vertex as V ^ v (p, q) = 
T)^ 1 + rJ^ M (p, q), with the SM and the beyond SM components given by 

= -gM v g, y (1) 
Tf u M (p,q) = ^[X(p-qg,u-p,q,)+ A'e^ (2) 

where A and A' are effective coupling strengths for, respectively, higher di- 
mension CP-even and CP-odd operators and we will assume that they are 
the same for W and Z bosons. These operators may be generated within the 
SM at higher orders of perturbation theory, although the resulting couplings 
are likely to be very small. In general, A and A' can be treated as momentum 
dependent form factors that may also be complex valued. However, we take 
the approach that beyond the SM (BSM) vertices can be generated from an 
effective Lagrangian, which treats A and A' as coupling constants The 



most striking difference between the SM and BSM vertices of eqs. (CQ) and 
(T51) is that the latter has an explicit dependence on the momentum of the 
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gauge bosons. It is this feature that is the source of the differences that the 
BSM vertices generate in the kinematic distributions of tagging jets in the 
VBF and HV processes, compared to the SM case. 

In the case of a spin-2 resonance coupling to VV states, we will follow 
Ref. 13J and adopt the following effective Lagrangian 

£ cff = \^ {f l9liV B^B^ + f 2 g, u W» a W^) (3) 

where T aj3 is the spin-2 field and B^ u , are the U(l) and SU(2) field 
strengths, respectively and A corresponds to a cut-off scale, which should 
be set at a value of the order of the TeV scale. The two operators in the 
equation above do not account for all possible terms of the most general 
spin-2-VV vertex and one could, for instance, include operators such as 
d ll d u T a p(Bt ia B l 'P + W^ a W u P) but which are of higher dimension. We re- 
strict ourselves to the terms proportional to fi and fi of eq. ()3]), setting the 
cut-off scale to A = 3 TeV in our numerical illustration . 

The purpose of our analysis is to identify observables and model inde- 
pendent tendencies in the kinematics of partons due to the anomalous HVV 
couplings above in electroweak processes in which the Higgs-like boson is 
produced in association with two jets, i.e. in vector boson fusion and in as- 
sociated Higgs production with vector bosons. For a given choice of the spin 
and parity assignments, we identify regions of the phase-space that would 
be populated differently by the Higgs-like boson with anomalous couplings 
compared to the expectation in the SM. 

In our analysis, the vertices for the Lagrangians in the SM and in BSM 
with spin-0 and spin-2 bosons are calculated in FEYNRULES [14] and passed 



to the program MADGRAPH 15], which is used for the generation of the matrix 
elements for Higgs production in VBF and VH. To obtain the cross sections 
and distributions at the hadronic level, the CTEQ6L1 parton distribution 



functions are used |T7j]. F° r completeness, associated SM Higgs production 
with two high pt jets in the gluon-gluon fusion process (ggF) is also consid- 
ered; the cross sections and kinematical distributions are obtained using the 



MCFM 16] program that incorporates the QCD process pp — > H + 2j + X at 
next-to- leading order (NLO). The factorization and normalization scales are 
set on an event-by-event basis to the transverse energy of the Higgs boson. 
For the selection cuts, partons are required to have transverse momentum 
p T > 20 GeV, rapidity \y\ < 20 and be separated by AR > 0.7. 
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3. Spin and parity determination in VBF 

The differences in the tensor structure of the different terms in the HVV 
vertex would significantly impact the kinematic distributions of the tagging 
jets produced in VBF. For example, in Ref. [lo| the difference in azimuthal 
angle between the scattered quarks (A^-) was found to be an effective dis- 
criminant for the various terms in eq. fl2]) for a scalar boson. As we will show 
below, the momentum dependence of the BSM vertex has further strong 
implications for a number of hadronic observables. These include the trans- 
verse momenta of the tagging partons^ {pTjiiPTj.2) and the rapidity difference 
{Ayjj) between them. The obvious advantage of looking at such observables 
is that they do not require a full reconstruction of the Higgs-like boson from 
its decay products and one can use all the search channels at our disposal at 
the LHC, including the H — > 77 and H — > t + t~ modes. 

In this section we describe the effect on the rapidity and transverse mo- 
mentum distributions of the jets due to each of the operators of eqs. (2-3), 
without considering any interference between them. In doing so it is impor- 
tant to understand that the distributions are independent of the strength or 
sign of the couplings (A, A', jfij/b) which only affect the value of the total 
cross-section. The strength of these couplings play an important role when 
one considers the simultaneous presence of these operators in the HVV ver- 
tex and we will discuss this at the end of this section. 

Figured] displays the difference in rapidity between tagging partons (Ay^-) 
for each of the higher dimensional operators. We see that the rapidity differ- 
ences for each of the operators are shifted to smaller values. While the BSM 
+ and 0~ cases display a behavior in this observable that is almost identi- 
cal, for a spin-2 particle the peak in the Ai/jj distribution is shifted to even 
smaller values. The dramatic behavior of this observable makes it a strong 
candidate to distinguish between the SM and the BSM vertex structures. 

We also observe that the transverse momenta of the tagging partons be- 
come significantly larger for pure BSM + and 0~ states. This feature is 
more pronounced for a spin-2 state. However, the Mjj distributions are not 
altered and we will thus not display them. 

The above mentioned features of the Ayjj and px distributions are due to 
the presence of momentum dependent structures in the vertices of eqs. (|2|3|) . 



1 Tagging partons are defined as the scattered quarks in the process pp — > Hjj. At 
NLO these are defined as the leading partons in transverse momentum. 
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Figure 1 : Normalized distribution of the difference in rapidity between the scattered quarks 
in the VBF process for each of the SM and BSM operators individually. 

The factors involving four momenta in the BSM vertices mean that the mo- 
mentum flow through the vertex allows for a greater push being given to the 
jets resulting in a shift of the vt distribution to larger values. This effect has 



been observed in Refs. |18|, ll3[ for both the scalar and the spin-2 case. 

An analytic calculation of the matrix element reveals that for both the 
SM and BSM cases, the vector boson propagator can be approximated to 
PTji PTj2^~ Ayji ) 2 for large values of the incoming parton momenta. The 
effect of this term is to push pj-j to smaller values and, at the same time, 
push Ai/jj to larger values. For the SM, the remainder of the terms in the 
square of the matrix element are proportional to s prji Prj2 Cosh(Ayjj), 
which further reinforces the occurrence of large rapidity differences. 

For of eq. ([2]), the additional terms have a dependence of the form 
(s pxji PTj2) 2 - This leads to much larger and flatter p T j distributions as 
compared to the SM. Although the rapidity dependence for these additional 
terms is complicated, it can be shown to have an opposite behavior to the 
propagator terms, unlike the SM case, pushing the rapidity difference to 
smaller values. Thus, the correlations between Aj^- and prju PT32 and/or 
Mjj critically depends on the tensor structure of the HVV vertex. 

This is further illustrated in Fig. [21 where we use the observable y/prji ■ prj2 
in view by the above considerations. This figure also suggests a new region 
of the phase-space for the exploration of physics beyond the SM in the Higgs 
boson sector by using hadronic observables. BSM vertices tend to populate 
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Figure 2: Two-dimensional distributions of the difference in rapidity and the square root of 
the product of the transverse momenta of the tagging partons. The upper plots correspond 
to the SM with electroweak (left) and ggF+2j (right) production while the lower plots to 
an admixture with a BSM + state with A = 1 (left) and a pure 2 + state (right). 

a region of the phase-space where the production of SM Higgs boson is sig- 
nificantly depleted. Three distinct regions of the phase-space in the Higgs 
boson production in association with two highest px partons are identified: 

• The QCD region in which the production of ggF+2j is dominant. Here, 
the partons tend to be close in rapidity and display moderate px with 
< JVTjx • VT32 >~ 50 GeV. In this corner of the phase-space, associated 
Higgs production with a vector boson with V — > jj contributes as a 
sub-leading process, for which < Jpxji • VTji >~ 45 GeV. 
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SM EW Hjj 


BSM 0+ 


o- 


2+ 


SM ggF+2j 


Process 
Acceptance 
a (fb) 


VBF VH 
0.06 0.04 
0.14 0.04 


VBF VH 
0.59 0.12 
1.43 0.29 


VBF VH 
0.55 0.18 
1.35 0.43 


VBF VH 
0.93 0.75 
2.25 0.86 


0.10 
0.27 


„h!W l„EW 
"BSMl"SM 




8.8 


8.7 


17.4 





Table 1: Acceptance, expected cross-sections (in fb) and the ratio of BSM Higgs boson 
cross-section to that in the SM populating the BSM electroweak region (Ayjj < 4 and 
Jprji ■ PTj2 > 100 GeV). Pure BSM states are considered. The cross-sections are given 
for the H — > 77 decay and 8TeV center of mass energy (see text). 



• The SM VBF region which is populated by the SM Higgs boson pro- 
duced via the VBF process. This corresponds to intermediate px tag- 
ging partons with < Jprji • Vrji >~ 55 GeV with a large separation 
in rapidity of < Ay 3 j >~ 4.5. Production of additional central partons 
is depleted, which is a typical signature of electroweak processes. 

• The BSM electroweak region which populates a corner of the phase- 
space defined b}@ ■Jprji • VT32 > 100 GeV and Ayjj < 4. The VBF and 
VH productions can be separated by requiring that the invariant mass 
of the leading partons be away from the weak boson invariant mass. 
A distinct BSM VBF region is defined where electroweak-like gluon 
radiation is expected to occur. This feature is important to further 
differentiate and suppress QCD Higgs boson production. 

Table [Q displays the acceptance of the SM Higgs boson and other BSM 
pure states in the BSM electroweak region. The latter is defined as Ayjj < 4 
and y/PTji ■ PTj2 > 100 GeV, where tagging partons are required to have 
Pt > 25 GeV and \y\ < 5. The table shows the acceptance of each pro- 
cess, separating the VBF and VH production mechanisms. The BSM cross- 
sections are assumed to be the same as the one for electroweak SM Hjj 



production which is normalized to the NNLO-QCD value [191 ]; all cross sec- 
tions are given for the H — Y 77 decay 

The SM Higgs boson production displays acceptance^] of less than 10%. 



2 The separation in rapidity in the BSM electroweak region tends to be somewhat larger 
compared to that displayed by ggF+2j production. For instance, < Ayjj >— 2.1 for a 
pure BSM + state compared to < Ayjj >— 1.7 for ggF+2j. 

3 The acceptance of the ggF process is reported with respect to the generator cuts given 
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In this process, additional gluons are radiated more copiously than in the 
electroweak processes. The fraction of the events in the BSM electroweak 
region that survive a veto on additional gluons with pt > 25 GeV and \y\ < 5 
is 50~^25%. The last row gives the ratio of the BSM contribution to that of 
the SM, where only electroweak processes are taken into account. In order 
to give a rough estimate of the signal significance for different scenarios with 
the if 77 decay, the non-resonant production of two photons and two 



partons is modeled with Madgraph [15j. The cross-sections of pure BSM 



signals, of background and the corresponding signal significance are scanned 
function of the lower bound on ^fprji • Prj2- 

It is found that for pure BSM scalar states maximal significance is ob- 
tained with y/pxji ■ VTj2 ^ 100 GeV and for a pure spin-2 state is obtained for 
y/PTji ■ PTj2 > 300 GeV. Assuming 8 TeV center of mass energy and 25 fb _1 
of integrated luminosity, a significance of approximately 1.4 a and 5 a can 
be achieved for the pure BSM scalar and spin-2 cases considered here, re- 
spectively. This sensitivity is achieved with signal-to-background ratios of 
0.08 and 1.5 for the pure BSM scalar and spin-2 states, respectively. This is 
defined in the di-photon invariant mass window 120 < M 77 < 130 GeV. 

It is important to note that the experiments at the LHC currently do not 
explore the BSM electroweak region. This region is vetoed by requiring a 
large separation in rapidity of the tagging partons, in order to isolate the SM 
VBF process. 

Another distinguishing feature is the relative increase of the cross-section 
with the energy of the collisions. This feature can be studied when the LHC 
begins its higher energy run. Figure |3] shows how the ratio of the cross- 
sections at -\/s = 8 TeV and yj~s = 13 TeV for the higher dimension operators 
differ from the SM case as the lower cutoff on the value of the di-parton 
invariant mass (Mjj) is changed (x-axis of the plot). The advantage of using 
these ratios, besides that they are free from many systematics [20j . is that 
they are independent of the values of the form factors A, A', f\ and / 2 . This 
ratio is greater than unity and rises with the Mjj cut for all the SM and BSM 
cases under consideration. For the spin-0 case this ratio is larger for the BSM 
operators as compared to SM. This occurs in spite of the total cross-section 
(without the application of any cuts) rising at the same rate (an increase by 
a factor of 2.3 from 8 TeV to 13 TeV) with yfs for I™ as well as each of 



in section. [5J 



9 




Figure 3: Left: ratio of the cross-sections of the VBF process at y/s — 13TeV to that 
obtained at y/s = 8 TeV as a function of the lower bound on the invariant mass of the 
scattered quarks Mjj ; results are shown for SM (black) and compared to those of the pure 
BSM + (red) and CP (blue) states. Right: the distributions of the rapidity difference of 
the scattered quarks for three different values of the upper bound on Mjj . 

the terms in T^ M . This result is a mere consequence of the difference in the 
energy dependence of the acceptance for different spin-parity combinations. 

As noted above for spin-2, the increase is much larger 7) and in fact the 
ratios of cross-sections grow much more rapidly for large values of the Mjj cut. 
This raises a concern about possible violation of unitarity for large Mjj values. 
However, there are several ways in which the spin-2 model can be unitarised. 
If the unitarisation prescription were to cut off the phase space concerned 
with large transverse momentum of the jets, then the Ai/jj distributions 
would shift to larger values simply due to the kinematic structure of the 
VBF mechanism. However, in a reasonable model, the unitarisation may 
also be implemented such that large values of Mjj are cut off. In this case, 
the peak in the At/jj distributions at small values will remain as shown in 
the right panel of Fig. |3j 

So far, our discussion was independent of the strength of the couplings in 
the HVV vertex. While this coupling strength is fixed in the SM, it is not 
the case for the BSM couplings and clearly this will affect the actual values of 
the cross-sections. Setting A = 0.4, A' = 0.47 and fx = 0.37 = f 2 ensures that 

a SM „ ^BSM^ 

An additional complication is worth noting. The acceptance 
for these operators are not the same in the same phase space, and one must be 
careful in correlating the number of events observed to a total cross-section. 
For example, applying an Mjj cut removes more events from the additional 
operators than from SM. One must be more careful with the total rate when 
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Figure 4: Normalized distributions of the difference in rapidity between the scattered 
quarks in VBF production: the SM Higgs case (solid lines) is compared with the case of 
BSM admixtures with a + term (left) and a 0~ term (right) as discussed in the text. 

looking at admixtures of these operators, which is possible only for the spin- 
case. The value of <t bsm may increase or decrease depending on the sign, 
magnitude and on whether the couplings are real or complex valued. 

Let us now consider the effect on the kinematic distributions, of the si- 
multaneous presence of the SM/BSM operators, which can happen only for 
spin-0. In Fig. HI the plot on the left compares the expectation from the SM 
(solid line) with an admixture with a BSM + scalar with A = 1 (dashed 
line) and A = i (dotted line). The plot on the right shows the same for the 
0~ scalar with A' = 0.1 -i and A' = 0.1. Only the strength of the coupling, 
and not the choice of BSM operators, affect this distribution. Hence, the 
plot on the right would look the same if a 0~ operator instead of the + one 
were used along with SM. Ay^ can therefore effectively be used to differ- 
entiate a SM from a momentum dependent HVV vertex, for larger values 
of these couplings. The plot on the left illustrates the fact that for not so 
small values of A/A', it is no longer possible to disentangle the SM Higgs bo- 
son from admixtures with BSM states using rapidity distributions. However, 
the angular correlation still shows small deviations from the SM case 

for A' = 0.1. One should note that the difference in the A<pjj distributions 
can be enhanced by looking at events with large pxj > 80 GeV. A similar 
observation is valid for comparable values of A. 

An important and interesting observation when considering admixtures 
of the BSM operators with SM is that for a certain choice of the parameters, 
SM = 1 and A = A' = 1, the A<j>jj distribution becomes indistinguishable 
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Figure 5: Invariant mass of the Higgs and the weak boson (left), the transverse momentum 
of the boson (center) and the separation in rapidity between the bosons (right). The 
solid line corresponds to the SM Higgs boson. The dashed, doted and dotted-dashed 
lines correspond to admixtures of the SM Higgs boson with a scalar CP-even state with 
A = 1,0.5 and 0.1, respectively. 

from the SM case. Although this is a special case and the cross-section is 
much larger here, it accentuates the importance of the Ayjj distribution and 
makes it necessary to complement them with A<pjj distributions. 

4. Kinematics in the VH mechanism 

Recently, associated Higgs production with a vector boson has been stud- 
ied to distinguish the case of a SM Higgs boson from pure BSM 0~ and 2 + 



states Ref. |21|. We survey in some detail effect of an admixtures of the SM 
operator with the BSM CP even and CP odd operator in case of VH pro- 
duction as well. We find that the invariant mass of the HV system provides 
a very good probe of simultaneous presence of SM and BSM operators. 

Figure [5] shows the invariant mass of the di-boson system, as well as 
the transverse momentum and the separation in rapidity between the Higgs 
and vector bosons. It compares the prediction for the SM Higgs boson and 
different admixtures with the CP even operator. With the enhancement 
of the boson transverse momentum, the overlap with the BSM electroweak 
region increases, as illustrated in Tab. [TJ 

5. Conclusions 

Probing the spin and CP structure of the recently discovered boson at the 
LHC is of prime importance. In this paper, the sensitivity to physics beyond 
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the SM via anomalous Higgs couplings to weak bosons is considered, where 
contributions from admixtures of BSM spin-0 and spin-2 states to the HVV 
coupling are tested. A phenomenological survey is performed by exploring 
various observables in the electroweak production of the Higgs boson with 
two high pt jets. It is found that the kinematics of the tagging jets in the 
vector boson fusion mechanism is strongly modified in the presence of mo- 
mentum dependent anomalous HVV couplings. In particular, compared to 
the SM, the correlation between the separation in rapidity and the transverse 
momenta of the scattered quarks varies significantly in the presence of new 
physics. The separation in rapidity is reduced while the transverse momenta 
becomes significantly larger; in addition, there is some complementarity be- 
tween the behavior of the azimuthal and rapidity separation between the 
scattered quarks. This defines a new corner of the phase-space that has not 
been explored by the ATLAS and CMS experiments to test the spin and CP 
structure of the Higgs to vector boson coupling. 

When this analysis was in the completion stage, Ref. [lij appeared on 
the archives, wherein the use of the rapidity separation Ay,j to distinguish 
between various CP and spin states is also advocated. In our paper, we have 
in addition investigated the effect of anomalous couplings on the correlations 
between y/PTjiPTj2 and Ayjj, thereby uncovering new regions of phase-space 
for exploiting the VBF production process. Besides showing that the ac- 
ceptance of the VBF process to the Mjj cut is different for different spin 
and parity assignments, we have evaluated how the acceptance (and hence 
the various kinematical variables) varies with the center of mass energy of 
the proton-proton collision and explored the use of relative rates between 
13TeV and 8TeV LHC to sharpen up the differences. Finally, in addition to 
the VBF process, we have also considered associated Higgs production with 
vector bosons, and have shown that some kinematical variables such as the 
invariant mass of the VH system are sensitive to anomalous couplings too. 
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